Thermodynamics of driven collisionless systems 
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A statistical theory is presented which allows to calculate the stationary state achieved by a driven 
system after a process of collisionless relaxation. The theory is applied to study an electron beam 
driven by an external electric field. The Vlasov equation with appropriate boundary conditions is 
solved analytically and compared with the molecular dynamics simulation. A perfect agreement is 
found between the theory and the simulations. The full current- voltage phase diagram is constructed. 



I. INTRODUCTION 

Unlike the equilibrium thermodynamics and statisti- 
cal mechanics, which are well developed after the pio- 
neering works of Boltzmann and Gibbs, our understand- 
ing of non-equilibrium thermodynamics is restricted to 
some special models and cases. Stochastic lattice gases 
have provided a fertile testing ground for studying non- 
equilibrium stationary states in driven systems [J, 0, 0| . 
These models exhibit a variety of phase transition aris- 
ing from a diffusive (collisional) relaxation. For some of 
these models local equilibrium and hydrodynamic equa- 
tions have been derived rigorously 

There are, however, other physical systems for which 
the approach to final stationary state is through a process 
of collisionless relaxation [j| [y, 0, 0, 0, OH Qj| . Gravita- 
tional systems and confined one component plasmas are 
just two such examples. For these systems the collision 
duration time diverges and the relaxation is governed by 
the collisionless Boltzmann (Vlasov) equation [ijj • In the 
thermodynamic limit, the collisionless relaxation process 
leads to non-Maxwell-Boltzmann velocity distributions, 
even for stationary states without macroscopic currents. 
Unlike normal thermodynamic equilibrium, the station- 
ary state which follows the collisionless relaxation de- 
pends explicitly on the initial distribution of particle po- 
sitions and velocities. In spite of this complication, it was 
recently shown that it is possible to construct a statistical 
theory that quantitatively describes these states [l(| EH • 

Beams of electrons driven by accelerating vacuum de- 
vices, like the thermionic valves, diodes, and magnetrons, 
also do not relax to the Maxwell-Boltzmann distribu- 
tion [l3| . Unlike the driven stochastic lattice gases, these 
systems, however, are intrinsically collisionless. An im- 
portant practical question concerns the kinetic tempera- 
ture distribution in thermionic devices in which the di- 
rected velocity produced by the electric field is compara- 
ble to the thermal velocity [l4j. This is particularly the 
case for the transitional region between Child-Langmuir 
and no-cutoff regimes in magnetrons, where the electric 



potential becomes comparable to the thermal energy [Tj] . 
Even when the final directed velocity is larger than the 
thermal velocity, there is a region near the emitting cath- 
ode where thermal effects are important. It is of great 
practical interest to determine the extent of these regions 
[T5L f]~6l| . Furthermore, since in these systems the collision 
duration time diverges, there is no local equilibrium, and 
one can not a priori postulate an equation of state re- 
lating the beam density and the beam temperature, as 
for adiabatic or isothermal processes. Instead, given the 
properties of thermionic filaments — such as say the ve- 
locity distribution of the emitted electrons — one should 
solve the boundary value problem posed by the Vlasov 
equation. The purpose of this Letter is to develop a theo- 
retical framework which will allow us to study the relax- 
ation dynamics and the stationary states of collisionless 
driven systems. 

As a prototype of a collisionless driven system, we con- 
sider a beam of electrons, accelerated by an external elec- 
tric field, traveling from an emitting (planar) cathode to 
a collecting (planar) anode across the device gap. The 
cathode, located at position x — 0, is kept at electrostatic 
potential ip(x — 0) = and is heated to temperature T c , 
resulting in the emission of electrons. After traversing 
the device gap, these electrons are collected at the cold 
anode (T a w 0) located at x = L and kept at potential 
(p(x = L) = V > 0. During the steady state operation, 
the region between the cathode and anode contains a to- 
tal of N electrons, resulting in a current density j. Our 
goal is to relate j, to the potential difference V, the num- 
ber of electrons N, the device width L, and the cathode 
temperature T c . 

For planar electrodes, particle distribution transverse 
to the x axis can be taken to be uniform. Further- 
more, the one particle distribution function for a colli- 
sionless system in a steady state must satisfy the sta- 
tionary Vlasov equation, 
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where e is the elementary charge, to is the electron mass, 
and / = f(x, v) is the static distribution function. In the 
thermodynamic limit, Vlasov equation becomes exact for 
particles interacting by long range potentials [17] . 

It can be readily seen that the distribution functions of 
the form f(x, v) = f[s(x, v)], where e is the mean particle 
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energy, e = mv 2 /2 — etp(x), satisfy Eq. ([T]). Therefore, 
if /(e) is specified at x — 0, / is then also determined 
for any other position, provided that the electrostatic 
potential <p(x) is known. This potential can, in turn, 
be calculated self-consistently from the solution of the 
Poisson equation 



d 2 <p(x) 
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where the particle density n(x) is given by n(x) = 
J v f(e)dv = J v f(mv 2 /2 — e<p(x))dv, the total particle 

number is N — A J v J Q f dxdv, and the transverse cross 
sectional area of the essentially ID device is A. To rep- 
resent both the thermal distribution near the cathode, 



and the fact that only particles with positive velocities 
actually move into the device gap, we choose at x = a 
unidirectional Maxwellian distribution of the form 
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where ks is the Boltzmann constant, T c is the cathode 
temperature, and no is the beam density at the cathode 
after the stationary state is achieved. The value of no can 
only be obtained once the full problem has been resolved. 
The distribution function over the length of the whole 
diode is then 
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Integrating the distribution function f[e(x,v)] over 
the possible values of velocity, we arrive at a nonlinear 
integro-differential equation for the electrostatic poten- 
tial, 
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dx 



It is important to note the difference between this equa- 
tion and the Poisson-Boltzmann equation obtained for 
usual collisional plasmas and electrolytes in the mean- 
field limit [l8j]. Equation §5§ can be solved numerically, 
to yield the electrostatic potential and the distribution 
function for the electron beam in the stationary state. 

For systems with long range interactions, Vlasov equa- 
tion should become exact in the thermodynamic limit. 
To confirm this for our system, we have performed molec- 
ular dynamics simulation of an equivalent one dimen- 
sional model. The simulated system consists of N s mu- 



tually interacting charged sheets of area A — each con- 
taining n s electrons of the same velocity — moving along 
the x axis, under the action of the external electric 
field produced by the grounded cathode ip(0) = and 
an anode kept at a fixed potential f(L) = V. The 
interaction potential between the two sheets G(xi , Xj ) 
is the Green's function of the Laplace equation, 
d 2 G(x,y)/dx 2 = 1/L S(x — y) with the boundary con- 
ditions G(x = 0,y) = G(x — L,y) — 0. Solving this 
equation we obtain G(xi,Xj) = x < /L (x > /L — l), where 
x < and £> are the smaller and the larger of the two par- 
ticle coordinates Xi and Xj. The effective Hamiltonian 
for the sheet dynamics is then 
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where e s = n s e and m s — n s m are the charge and mass 
of each sheet respectively. The acceleration of each sim- 
ulated sheet then follows from the canonical equations of 
motion, 
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where 7ii /e ^*(" s,l *'l is the number of sheets to the 
left (right) of the one considered, and x denotes the posi- 



tional average x = J^j X j/N s - Since < x < L, from Eq. 
d?]) one sees that the electron acceleration at the device 
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entrance where n\ le ^^ — > and n' 1 " 9 ' 1 ^ — > iV s satisfies 
eV/mL — Ne 2 /2e mA < v{x = 0) < eV/mL, which 
reveals that in space- charge dominated devices where 
eV/mL < Ne 2 /2eomA, acceleration at beam entrance 
may be zero or even negative [2(| • When the acceler- 
ation vanishes, the associated current is denoted as the 
limiting one. Since we wish to describe a hot cathode 
and a cycling current inside the device, we adopt the 
following strategy. We advance the simulation in small 
time steps, always obeying Eq. ([7]). Whenever a particle 
crosses the anode and exits the system, it is re-injected 
at the cathode position. At this point all the particles in 
a small region S s around the cathode are re-thermalized, 
so as to ensure that the distribution there keeps its orig- 
inal form of a truncated Maxwellian. The width S s must 
be sufficiently small, S s <C L, but apart from this condi- 
tion its precise value is arbitrary. The simulations were 
performed with 5 S /L = 0.01 and N s = 50,000. In all 
cases we start with a uniform distribution of sheets and 
compute the observables only after the system reaches its 
final stationary state. 

To compare the predictions of the theory with the re- 
sults of the simulations, we consider the density and the 
temperature distributions inside the diode. The kinetic 
temperature is defined as 

k B T{x) =^{x) ~v 2 {x). (8) 

where the over-bar denotes the velocity average at a given 
position x. The theoretical averages are calculated using 
the distribution function f[s(x, v)], while in the simula- 
tions, the averages are performed over the particle ve- 
locities within narrow bins along the x axis. Note that 
because of the asymmetry of the velocity distribution at 
x = 0, T(0) + T c . 

It is convenient to scale space and time with the diode 
length L and the plasma frequency lo 2 = Ne 2 /cqitiLA, 
respectively. Dimensionless coordinate and velocity can 
then be defined as x* — x/L and v* = v/Lu p . In ad- 
dition, Eqs. ([7]) and ^ show that adimensional volt- 
age and adimensional temperature can be defined as 
V* = eV/mL 2 ui 2 and T* = k B T/mL 2 uj 2 , respectively, 
and serve as the control parameters for the system. 

In Fig. []Ja) the scaled temperature T* is plotted 
against the scaled coordinate x* . We consider T* = 0.05 
and also consider a device operating at its limiting cur- 
rent, ii(x = 0) = 0. A striking feature of this plot is that 
the temperature drops rapidly as one moves away from 
cathode towards anode. We next study the dependence 
of scaled density n*(x*) = n(x)AL/N along the length 
of the diode. The density is very high near the cathode, 
where the average velocity is small. It then drops rapidly 
towards the anode, where particles are accelerated up to 
high speeds, see Fig. [ljb). Agreement between the sim- 
ulations and the theory for both the kinetic temperature 
and density is excellent. 

We now study the current-voltage phase diagram of 
the device. In general, current is a function of the volt- 
age drop, the temperature, the gap length, and the total 
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FIG. 1: (a) Temperature and (b) density versus position in 
the case of the limiting current and T* = 0.05. Solid lines rep- 
resent the theoretical results, while the circles are the results 
of the simulations. 

charge of the device. However, by measuring the time 
in units of one over the plasma frequency u>~ , and the 
length in units of the gap length L, we can scale away two 
of these variables. The current density can be calculated 
using 

j = -e J vfdv. (9) 

Since in the steady state the current does not depend on 
either time or coordinate, integration along the x axis 
and over the cross sectional area A yields 

jLA=-e J vfdvdxd 2 r ± . (10) 

Furthermore, since the current density is measured in 
units [j] ~ [Nev/AL], rescaling it in terms of the gap 
length and the plasma frequency, we can write the re- 
duced current density as j* = ~ej / e^mLuj^, which then 
satisfies 

f=W{V*,T*), (11) 

where v* is the reduced velocity averaged over all the 
particles. The reduced average velocity, in turn, must be 
a function of the two previously introduced control pa- 
rameters: the reduced voltage and temperature. Eq. (|ll|) 
is in fact a similarity transformation relating systems 
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with different charge, length, temperature, and poten- 
tial difference. In Fig. [2] we plot j* vs. V* for vari- 
ous T* . The phase diagram provides all the information 
about the current-voltage characteristics for all possible 
planar diodes. The first feature to note is that all the 
different curves emanate from the limiting current back- 
bone, which traces a temperature dependent path in the 
j* x V* plane. To the left of the limiting current bor- 
der, indicated by the solid line in Fig. [2 the distribution 
function can no longer be described by a unidirectional 
Maxwcllian, such as the expression ^j. The transition 
resembles Bose-Einstein condensation (BEC). In the case 
of BEC — below the critical temperature — a macroscop- 
ically populated ground state appears, and only a frac- 
tion of particles remains in the excited states. Similarly, 
in the case of our diode, to the left of the limiting curve, 
part of the charge must be expelled from the system be- 
fore a stationary state can be achieved. As T* — > oo, the 
voltage effects become negligible compared to the ther- 
mal ones, and the beam density becomes uniform across 
the gap. In this limit, it is possible to show that the 
dimensionless backbone curve asymptotes to a vertical 
line, V* = 0.5, Fig. H 
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FIG. 2: Characteristic curves of j* vs. V* . The thick solid 
line represents the limiting current and the thick dashed line, 
the zero temperature limit. To the left of the solid curve, 
charge must be expelled from the system before a stationary 
state can be achieved. Dotted lines represent the theoretical 
results for the indicated temperatures. The circles are the 
results of the simulations at the same temperatures. 

To conclude, we have studied the dynamics of collision- 
less driven systems. Unlike the stochastic lattice gasses 
which are significantly abstracted from reality, the mod- 
els studied in this paper are very similar to real electronic 
devises, such as the thermionic valves, diodes, and mag- 
netrons. Furthermore, differently from the lattice gases 
whose dynamics is diffusive, the distribution function of 
collisionless systems satisfies the Vlasov equation. For 
the class of driven systems introduced in this letter, the 
stationary state Vlasov equation can be solved exactly. 
The theory developed in this paper should, therefore, be 
relevant to the design and operation of real electronic 
devises. 

It is important to stress that in the absence of colli- 
sions, a charged beam does not relax to an equilibrium 
with a known equation of state. In fact, the thermody- 
namic temperature is defined only in the vicinity of the 
hot emitting cathode. Away from the cathode, dynamics 
is controlled by the collisionless Vlasov equation, which 
has to be solved as a boundary value problem. Once 
the solution is obtained, all the macroscopic quantities 
can be determined via appropriate averages. The kinetic 
temperature is found to vary strongly across the device 
gap, precluding the use of conventional isothermal or adi- 
abatic assumptions and of the hydrodynamic formalisms. 
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